function g = elliptic_g(x,y,caseNum,varargin)
global alpha theta;

        g = (x.^2+y.^2).^(1/alpha).*sin((2*atan2(y+eps,x+eps)+theta)/alpha);
   